Extension of relativistic dissipative hydrodynamics to third order 
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Following the procedure introduced by Israel and Stewart, we expand the entropy current up to 
the third order in the shear stress tensor Tv a/3 and derive a novel third-order evolution equation for 
ir°^. This equation is solved for the one-dimensional Bjorken boost-invariant expansion. The scaling 
solutions for various values of the shear viscosity to the entropy density ratio rj/s are shown to be 
in very good agreement with those obtained from kinetic transport calculations. For the pressure 
isotropy starting with 1 at tq — OAfm/c, the third-order corrections to Israel-Stewart theory are 
approximately 10% for rj/s = 0.2 and more than a factor of 2 for ij/s — 3. We also estimate all 
higher-order corrections to Israel-Stewart theory and demonstrate their importance in describing 
highly viscous matters. 



PACS numbers: 47.75. +f, 24.10. Lx, 24.10.Nz, 12.38.Mh, 25.75.-q 



o . 

<: 

(N 



43 

Of 
i ' 

or 

43 



(N 
> 
O 
O 

in 
o 

&\ 
O 



X 



A causal theory of relativistic dissipative hydrodynam- 
ics was first formulated by Israel and Stewart [l| and has 
been successfully applied to study and understand a wide 
range of phenomena observed in ultrarelativistic heavy- 
ion collisions The success may lie on the specula- 
tion that the viscosity of the hot dense matter created 
experimentally is really small. Since the Israel-Stewart 
theory is a second-order theory, which neglects the higher 
orders in viscous stress, it is natural to look for its limit of 
applicability. Especially for one-dimensional expansion 
with the Bjorken boost invariance Q the Israel-Stewart 
theory has led to a reheating of the expanding medium [6] 
and to a negative pressure pj], if the viscosity is large or 
equivalently the starting time of expansion is small due 
to the scaling behavior of the solutions. These effects are 
unphysical and thus show the breakdown of the Israel- 
Stewart theory. Moreover, in Refs. d, |{|, the scaling 
solutions of the Israel-Stewart hydrodynamic equations 
are compared with those obtained from kinetic transport 
calculations. The breakdown of the second-order the- 
ory has been found at rj/s > 0.5. From comparisons 
of dissipative hydrodynamics calculations to the elliptic 
flow measurements one infers that rj/s ~ 0.5 is an upper 
limit for the shear viscosity to entropy density ratio of a 
quark-gluon plasma (QGP) [TcHl^ ]. It is thus of interest 
to extend the second-order Israel-Stewart theory in order 
to include higher-order corrections and then to quantify 
their effects on observables. 

Using the Grad's 14-moment method and neglecting 
bulk pressure and heat flow, the off-equilibrium distribu- 
tion function is given 
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f(x,p) = / (x,p)(l 
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fo(x,p) is the Boltzmann distribution in kinetic equilib- 



fo(x,p) = (2?r) 3 Xexp 



T 
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where A denotes the fugacity describing the chemical 
equilibration, T is the local temperature and is flow 
velocity. The expression in Eq.([T]) is first order in 
= T^ v — which is the deviation of the energy- 
momentum tensor from its equilibrium form. We con- 
sider a system of massless particles which leads to a van- 
ishing bulk pressure. Using the Landau matching condi- 
tions u„u v T^ = u^N^ = we obtain C = 3/(8eT 2 ) 
Q. The Landau matching conditions are equivalent to 
e = e eq and n = n eq in the comoving frame, where 
e eq = 3\T 4 /ir 2 and n eq = XT 3 /tt 2 denote the equilib- 
rium values for energy and particle densities according 
to Eq.©. This allows us to define the effective temper- 
ature T = e/ (3n). 

We follow the approach introduced by Israel and Stew- 
art and use the entropy principle to derive an evolution 
equation for the shear tensor. The entropy current 
can be calculated according to the kinetic definition: 
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ln(/) will be expanded to the third order in 
Co-KnupPp" [see Eq.Q]. We obtain 
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where sq — — J d 3 p /n(ln/o — 1) = An — nlnA is the 
entropy density in kinetic equilibrium and fi% = 9/(4e). 
The value of ft% we obtained is just the same as derived 
in the Israel- Stewart theory [HQ. 

Up to second order in 7r Q/3 the entropy current in Eq. (TJ|) 
is the same as in Israel-Stewart's a pproachQ, 

How- 
ever, it looks different from the expression found in Refs. 
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[HI [l6j for conformal fluids. There, the entropy current 
is constructed to include all possible terms up to the sec- 
ond order in gradient of the flow velocity. In general, the 
shear tensor n a ^ can be expanded in terms of gradient 
of the flow velocity by introducing transport coefficients. 
Up to the second order the expansion is explicitly given 
in Ref. [13] for conformal fluids. If we insert this expan- 
sion into Eq.©, we cannot obtain the same expression 
as in Refs. [la, lltij . This probably indicates that our 
ansatz Eq.© could be modified. However, an extension 
of Grad's approach has not been developed so far. 

The main goal of this work is to investigate whether 
the third-order term cx 7r Q ,^7r"7r^' T affects local entropy 
density production. We expect that this higher-order 
term reduces the local entropy density s for typical off- 
equilibrium initial conditions in heavy-ion collisions, e.g. 
the color glass condensate, i.e. it acts in same way as 
the second-order term introduced earlier by Israel and 
Stewart. 

Taking the divergence of the entropy current and using 
the Gibbs-Duham relation [H, 0| , we obtain 
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side of Eq. ([7]) is always non-negative, whereas the sec- 
ond term is first negative (t^O > 1) and then positive 

M < !)• 

The entropy production Eq. ([5]) can be rewritten as 
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Note that each term in the brackets has to be traceless 
for further derivations and we have thus applied the pro- 
jector operator to the quadratic terms in 7r Q/3 . According 
to the second law of thermodynamics, the entropy pro- 
duction should be non-negative. This is guaranteed if we 
impose a linear relation between the dissipative flux 7r Q ^ 
and the expression in the square brackets, which can be 
interpreted as a thermodynamic force: 
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where 



(9) 



= V<"u"> - ( ^(A£A£ + A- Ap - ^A^A^J V 

(6) 

and A Q( g ~ g a p — u a up with the metric g a p — 
diag(l, — 1, — 1, — 1). ir al3 is the derivative with respect 
to r = y/t* - z 2 . 

First we demonstrate that the last term in Eq.© al- 
ways has a positive part, which is not relevant for the 
derivation of a relaxation equation for n a ^ . In the re- 
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laxation regime, which is determined by tt 
the sign of the product — 7r Q , ( g7r"7r' 3<T is opposite to that of 
— Tr a pTT"Tr l3 ' T . Thus the last term in Eq.© is positive. If 
the system is initially outside the relaxation regime, the 
sign of — TTapir^ir 13 ' 7 will be the same as — 7r a/ 37r"7r^ <T , i.e., 
negative. Then we can separate the last term of Eq. © 
into two terms: 
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where 6 — d^u^ denotes the inverse of the expansion 
scale and t„ denotes the intrinsic relaxation time scale 
on which the relaxation of the shear pressure toward the 
Navier-Stokes value sets in. t„9 is thus the local Knudsen 
number [8|, [l8(. We now argue that for t^Q > 1, the sys- 
tem is outside the relaxation regime; whereas for t v 8 < 1 
relaxation sets in. Thus the first term on the right-hand 



where ij is the shear viscosity. Note that if the sys- 
tem is initially in the relaxation regime the last term 
in Eq.© should not be included in the square brackets 
since it is a positive contribution to the entropy produc- 
tion [see Eq.©]. By separating Eq. © into two non- 
negative parts, we have assumed the maximum entropy 
production, which corresponds to the fact that interac- 
tions among matter constituents tend to drive the system 
toward equilibrium as fast as possible. Dividing both 
sides of Eq.© by 2r7/3 2 and using the Israel-Stewart defi- 
nition of the relaxation time 7v = 277/32 we finally obtain 
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Equation (|10[) presents a novel third-order evolution 
equation for the shear tensor. This is the main result 
of this work. We note that T n ~ tv/t is of the same 
order as ir a ° /T 4 for large r. Therefore, the last term in 
Eq. (fTT)]) is of fourth order in n a ^ /T 4 for large r and 
should be taken off from a certain time. More discussion 
on the higher-order corrections will be given later in this 
article. 

Neglecting the last two terms in Eq. (TTU1) . the second- 
order Israel-Stewart equation is recovered. As men- 
tioned before, it does not contain all second-order terms 
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found in recent works [!, [l6|, [T3, ■ Considering one- 
dimensional boost-invariant expansion, the only differ- 
ence is the second-order term ^-tt^tt 17 ^ [l7| . which 
does not appear in our equation. However, recently it 
has been shown in [?J that including this term does not 
help prevent unphysical behavior. We will demonstrate 
below that this behaviour will be corrected by the third- 
order term found by us, which has a form similar to 
■^tTt^tt' 7 ^ but contains an additional gradient. 

In the following we explicitly give the third-order 
viscous hydrodynamic equations for a one-dimensional 
boost-invariant expanding system. In this particular 
case the heat flow vanishes as we assumed at the be- 
ginning. The four flow velocity is u M = ^(t, 0,0, z) and 
thus 9 — <9 M w M = 1/t. In the comoving frame the shear 
tensor is diagonal, with a positive shear pressure 7r > 0: 



tt^ = diag(0, tt/2, tt/2, -tt). 



(11) 



The evolution equation for the energy density e and the 
number density n follow from the conservation of the 
energy momentum tensor, d v T v0 — 0, and the particle 
number conservation, d^N^ = 0: 
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Using Eqs. (112]) we obtain 
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according to Eq. (ITU1) in the comoving frame. 

To count the orders in n/e we multiply both sides of 
Eq. (Till)) by Tv/e and obtain 
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We realize that at late times when the relaxation toward 
equilibrium starts, t„/t, which is the local Knudsen num- 
ber, is of the same order as 7r/e, which can also be inter- 
preted as the inverse of the Reynolds number 6j. Thus, 
the last term in Eq. (TTJ} is of fourth order and has to 
be omitted for consistency. At early times tv/t > 1 all 
higher-order contributions have to be included. For sim- 
plicity we neglect the last term in Ea. (fT4)) for all times 
and obtain 
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The term — 37r 2 /(er), which gives the difference from the 
Israel-Stewart theory, additionally damps the increase of 
the shear pressure w. 

Comparing with Eq. (|15[) . one can assume that the 



higher-order terms are of the form - 
we take a heuristic expression 
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In practice, 



(16) 



where the number x is supposed to contain approximately 
all corrections of and beyond the third order. A simi- 
lar approach, an improved hydrodynamic theory, which 
includes higher-order gradient terms, has been recently 
proposed [201. Equation (Til))) should be valid even for a 
free streaming of particles, which is an extreme case with 
an infinite shear viscosity. Thus 7V ~ 77 = 00. For a 
one-dimensional free streaming with the Bjorken boost 
invariance the energy density decreases as e = — e/r, 
which leads to n = e/3 according to Eq. (fT2")l . The 
latter indicates the vanishing of the longitudinal pres- 
sure. Putting these into Eq. p^|) we obtain x = 5/3. 
We see that the damping of tt becomes weaker when all 
higher-order contributions are taken into account. This 
also indicates that the higher-order corrections may have 
an oscillating behavior. 

Using Eqs. (|12p and (fTBjl we obtain 
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For x — 0, which indicates the transition to the Israel- 
Stewart theory, the derivative of ir/e at ir = e/3 can 
be positive for sufficiently large t„. (or 77). This means 
that 7T can be larger than e/3, which leads to a negative 
longitudinal pressure and thus, is unphysical. Note 
that the negative effective pressure phenomenon has been 
observed in Ref. Q using the full second-order equations 
presented in Refs. [!, [l9j]. In contrast, for x — 3 or 
x = 5/3 in Ea.dTB)) the derivative of ir/e is negative at 
tt — e/3. We realize that the third-order (and higher- 
order) corrections prevent the unphysical behavior that 
appears in the second-order theory. 

To demonstrate the significance of higher-order correc- 
tions to the Israel-Stewart equation we solve the hydro- 
dynamic equations (Ti"2"j) and (|16l) with x = (IS), x = 3 
(third order), and x = 5/3 (all orders approximation), 
respectively. We also compare these solutions with those 
calculated from a transport model, the parton cascade 
B AMPS [2l[ , which has recently been applied to investi- 
gate a wide ran ge o f phenomena such as the buildup of 
the ellip tic flow [22| , the energy loss of high-energy glu- 
ons [23j, the extraction of the second-order viscosity co- 
efficient fej , and the formation and propagation of shock 
waves [18| in ultrarelativistic heavy-ion collisions. 

In this work, we implement only elastic collisions with 
an isotropic differential cross section. In this case, the 
shear viscosity coefficient can be calculated via [8[ 



6 T 
5 CT22 



(18) 



Using the equilibrium entropy density s = An — n In A we 
obtain the cross section as a function of 77/s: 
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We adjust the cross section locally to keep 77/s constant in 
the BAMPS calculations. In hydrodynamic calculations 
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FIG. 1: (Color online) Time evolution of the pressure isotropy 
for various rj/s values. Symbols present the results of BAMPS 
calculations. The solid, short dashed and long dashed curves 
show the solutions of the Israel- Stewart theory (x — 0), of 
the third-order hydrodynamics (x = 3), and of the heuristic 
consideration of all-orders contributions (x = 5/3). x is the 
parameter in Eq. (fl6)l . 

the shear viscosity is given via r\ = ( (An — n In A). The 
initial condition of the one-dimensional boost-invariant 
expansion is assumed to be a thermal state (A = 1) with 
an initial temperature of To = 500 MeV at the initial 
time To = 0.4 fm/c. 

Figure [1] presents the time evolutions of the pressure 
isotropy pl/pt = (p — n)/{p + 7r /2), where p = e/3, 
for various constant rj/s values. In all hydrodynamic 
and transport calculations, the pressure isotropy first de- 
creases from the equilibrium value 1 and then turns to 
increase toward the equilibrium value. The time scale of 
the minimum of the pressure isotropy is proportional to 
the relaxation time scale r n = 2?7/32, which is approxi- 
mately proportional to the rj/s ratio. 

Second, the larger the rj/s ratio, the wider the three 
hydrodynamic solutions go apart. Whereas for r\/ s = 
0.05 the curves are almost identical, for rj/s = 3 the third- 
order and approximative all-orders treatments bring a 
correction of 300% and 200%, respectively, to the result 
of the Israel-Stewart theory. [For rj/s = 1 (not shown) 
these corrections are 200% and 100%, respectively.] The 
higher-order corrections are essential for prevention of 



appearance of unphysical negative pressure, as observed 
in the IS solution for rj/s = 3. 

Third, comparing the results (symbols) from the 
BAMPS calculations with the solutions (solid lines) of 
the Israel-Stewart hydrodynamics, we see a perfect agree- 
ment for rj/s = 0.05 and a 10% difference for rj/s = 0.2. 
For rj/s = 0.4 the difference is already more pronounced, 
30%; and for r}/s = 3, it is significant - almost 200%. 
The second-order Israel-Stewart theory breaks down for 
large rj/s values. This is in line with the findings in Refs. 

Finally, the hydrodynamic solutions in the third-order 
and approximative all-orders treatments are much closer 
to the kinetic transport solutions than the second-order 
ones. The heuristic ansatz for including all orders cor- 
rections leads to remarkably good agreements with the 
BAMPS results for rj/s values in a wide range, although 
it is hard to extend this treatment to a general form such 
as Eq. (fTTJjl . We see that the main deviation of the third- 
order results from the transport ones appears during the 
early times, at which the local Knudsen number t^/t is 
larger than 1 and the dissipative effect is still strong, so 
that the third-order treatment is not sufficient to describe 
the real hydrodynamic expansion. This is also the rea- 
son why the approximative all-orders treatment provides 
a somewhat better description at early times. 

If we use a smaller initial time To, which will be the case 
at the CERN Large Hadron Collider (LHC), the system 
would evolve faster out of equilibrium. Smaller rj/s val- 
ues have to be chosen to obtain the same quantitative 
behavior as in Fig{T] 

In summary, we have derived a novel third-order evo- 
lution equation for the shear stress tensor in dissipa- 
tive hydrodynamics. The higher-order corrections to 
the second-order equation from the Israel-Stewart the- 
ory have been investigated. We found that the third- 
order contribution is essentially needed to prevent un- 
physical behaviors, which may occur for large rj/s ratios. 
Comparisons between hydrodynamic and kinetic trans- 
port calculations have quantitatively demonstrated the 
significance of higher-order corrections for description of 
expanding viscous matters. The solutions of the higher- 
order hydrodynamic equations showed good agreement 
with the kinetic transport results in a wide range of the 
■q/s ratio, although the derived equation does not contain 
all possible second-order terms obtained in the literature. 
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